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— *  A  simple  homogeneous  rock,  Arkansas  novaculite,  was  chosen  for  dynamic  tensile 
failure  studies.  Experimental  techniques  for  performing  dynamic  tensile  tests  on  rocki 
were  developed  and  used  to  provide  meaningful  measurements  of  dynamic  tensile  strength 
A  way  to  eliminate  the  occurrence  of  undesirable  radial  cracking  of  the  specimen  durinj 
testing  was  found;  specimen  recovery  and  fracture  analysis  methods  were  implemented. 

Novaculite  specimens  were  characterized  with  respect  to  static  and  dynamic  fractw  e 
strength,  fracture  toughness,  defect  and  microstructure,  and  other  properties  and 
parameters  pertinent  to  the  fracture  behavior^  )A  particularly  important  accomplish¬ 
ment  was  the  quantitative  determination  of  th^size  distribution  of  inherent  cracklike 
defects,  for  this  is  thought  to  play  a  dominant  role  in  the  comminution  characteristics 
of  the  rock,  and  is  thus  a  necessary  parameter  for  a  dynamic  fracture  model.  The  quasi 
Static  tensile  strength  was  determined  by  means  of  expanding  ring  tests  to  be  440  ± 

20  kg/cm  .  The  tensile  fracture  strength  of  govaculite  was  measured  in  high  loading 
rate  gas  gun  experiments  to  be  410  ±  20  kg/cin  and  found  to  be  independent  of  defect 
orientation.  Thjjs  th|  tensile  fracture  strength  is  strain  rate  insensitive  in  the 
range  10  to  10  sec  .  A  Griff ith-Irgin  tyye^relation  is  suggested  for  the  failure 
criterion  unger  dynamic  loads  <yf  =  10  c^  "  where  <y  is  the  fracture  strength  in 

kilograms/cm  and  c  is  the  maximum  cracic  radius  in  centimeters. 
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ABSTRACT 


This  report  presents  the  results  of  the  first  year  of  a  three- 
year  experimental  and  theoretical  program  to  study  dynamic  tensile 
failure  of  rock.  Techniques  for  loading  rock  specimens  in  tension 
at  high  rates  in  a  gas  gun,  and  recovering  them  are  described.  Ytterbium 
stress  gages  were  used  to  determine  loading  histories  and  measure  dynamic 
tensile  strengths. 

All  experiments  were  carried  out  on  Arkansas  novaculite,  a  quartzite 
chosen  for  its  simple  homogeneous  structure.  The  fracture  mechanism  is 
inferred  from  the  results  of  f ractographic  examinations  of  recovered 
specimens , and  a  computational  model  for  dynamic  tensile  fracture  in 
rock  based  on  the  hypothesized  fracture  mechanism  is  proposed  and 
partially  developed.  Future  work  is  aimed  at  further  development  and 
generalization  of  this  model  to  allow  predictions  of  fragment  size 
distributions  in  rock  from  known  dynamic  tensile  loading  histories. 
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‘  5  i 

INTRODUCTION  , ! 

i  *  \ 

Fracture  in  rock  under  high  rate  tensile. loading  conditions  is 
more  complex  and  considerably  less  well  understood  than  fracture  under  | 
quasi-static  conditions.  There  is  at  present  no  satisfactory  theoretical 
basis  for  predicting  dynamic  failure  behavior,  although  the  advantages  i 
of  having  such  a  basis  are  many.  An  understanding  of  rock  failure  under 
dynamic  tensile  loads  would  be  most  useful  in, the  solution  of  practical! 

i 

mining  and  civil  engineering  problems .  With  such  knowledge  rapid 
excavation  could  be  done  more  safely  and  economically,  the  stability  of. 
structures  in  rock  could  be  designed  and  evaluated  with , more  confidence, 
and  the  efficiency  of  rock  disintegration  processes  could  be'  improved. 

It  is  thus  the  objective  of  this  three  year  program  to  develop  a  model 
for  rock  fracture  which  can  be  used  to  predict  failure  behavior  under 

i 

high  rate  tensile  loading.  .  , 

This  report  describes  the  progress  made  during  the  fi,rst  year.  In 

this  section  we  discuss  the  strength  of  rockB,  compare  their  frapture 

! 

behavior  under  quasi-static  and  high  rate  loads,  and  review  previous 
attempts  to  model  dynamic  fracture  of  materials  in  general.  Properties1 
and  microstructural  features  of  Arkansas  novaculite,  the  specimen  material, 

i 

are  given  in  Section  II,  the  results  of  high  rate  loading  experiments  are 
presented  in  Section  III,  a  fracture  mechanism  is  proposed  in  Section >IV, 
and  a  dynamic  fracture  model  is  outlined  in  Section  V.  Plans  for  improving 
and  substantiating  this  model  are  stated  in  Section  VI. 

i 

The  Quasi-Static  Strength  of  Rocks  , 

Generally  speaking,  there  is  no  well-defined  yield  strength  or 

I 

fracture  strength  for  rock  materials,  Strength  values  measured  on  ten  or 

twenty  similar  rock  samples  may  vary  by  more  than  an  order  of  magnitude, 

i  i  . 

even  when  extreme  care  is  taken  to  ensure  that  the  samples  are  as  much 

1  I 

alike  as  possible.  Therefore,  rocks  cannot  in  general  be  characterized 


i 


I 


by  a  single  strength  Value, as  is  usually  possible  with  metals,  and  it 

'  ■  •  i 

is  of  little  help  to  use  a  criterion  such  as  "Fracture  occurs  when  the 
fracture  strength  of  the  material  is  exceeded.", 

1  i  1  1 

1  .1  f  .  r 

The  main  reason  for  the  large  scatter  m  the  measured  strengths  of 

*  ■  i  . 

rocks  is  most  probably  the  wide  range  in  the  sizes  of  inherent,  internal 

,  .  p  J 

cracks.  It  is  well  established  that  preexisting  cracks  control  the 
strength  6f  brittle  or  semibrittle  materials  and  that  the  larger  the 
cracky,  the  weaker  the  material.  In  particular,  the  size  of  the  largest 
crack  determines  the  <juasi-static  fracture,  strength.-  Rocks  typically 

have  a  very  pronoupced  defect  structure,  whose  size  range  is  much  broader 

! 

-than  that  typical  of  metals. 

Deliberate  processing  procedures  for  metals  control  the  maximum 

i  . 

flaw  size;  hence  the  fracture  strength  of  metal  is  rather  constant  from 

quantity  to  quantity.  fn  contrast,  the  haphazard  conditions  under  which 

1  i 

naturally  occurring  rock  is  formed  result  in  a  wide  distribution  of  flaw 
sizes  and  corresponding  scatter  in  fracture  strengths.  Thus  in  formulating 
a  predictive  capability  for  tensile  failure  in  rocks,  it  is  necessary  to 

take  into  account  the  size  distribution  of  inherent  cracks. 

1  2  1  1 

The  Griffith  theory  links  the  sample  strength  with  the  size  of  the 

maximum  crPck  and  provides  the  most  ireliable  expression  for  predicting 

,  l  • 

quasi-static  fracture  in  rock;  The  two  main  assumptions  of  the  theory, 

(1)  that  the  material-  is  linear  and  elastic  and  (2)  that  the  material 
contains  a  population  of  flaws,  are  well  satisfied  in  most  rocks. 

Griffith  related  the  brittle  fracture  stress  g  to:  the  elastic  modulus  E, 
the  surface  energy  y,  and  the  radi,us  of  the  largest  crack  c,  in  the 
foliowing  way: 


O’  i= 


* 

where  k  is  a  'geometrical  constant  of  'the  order  of  one. 

-  i  i 

*  k  does  in  fact  Contain; the  Poisson  ratio,  but  since  variations  in  this 
,  property  from  rock  to  rock  ar4  usually  insignificant,  k  may  be  considered 
to  depend  only  on  specimen  geometry. 

i  2 

•  I  ,  ! 


I 


Since  all  material  properties  are  contained  in  the  product  (Ey)^,  this 
term  characterizes  the  response  of  the  material  to  fracture.  Specifically, 
(Ey)^  is  the  material  property  that  indicates  the  basic  propensity  of  the 
material  for  brittle  tensile  failure.  In  fracture  mechanics  terminology^*^ 
it  corresponds  to  the  critical  stress  intensity  factor  K 

Ic 


<7 


*  K  /c 
Ic 


h 


This  inherent  fracture  resistance  may  be  considered  a  material  property, 

but  it  is  not  necessarily  a  material  constant.  In  general,  its  value  is 

a  function  of  geometry,  strain  rate  and  temperature. 

According  to  the  Griffith  theory  or  the  concepts  of  fracture  mechanics, 

the  quasi-static  strength  of  a  rock  sample  can  be  predicted  knowing  (1)  the 

inherent  fracture  resistance  of  the  rock  (Ey)^  or  K  and  (2)  the  size  of 

Ic 

the  largest  most  favorably  oriented  flaw  in  the  sample.  The  inherent 

fracture  resistance  may  be  evaluated  either  by  measuring  the  elastic  modulus 

and  the  fracture  surface  energy  individually  in  separate  experiments,  or 

by  measuring  the  strength  of  a  specimen  in  which  the  maximum  flaw  size  is 

known.  Statistical  methods  are  usually  employed  to  determine  the  maximum 

5,6 

flaw  size.  A  Scheil-type  statistical  transformation  of  flaw  traces  on 
sectioned  surfaces  will  be  shown  to  have  particular  application  to  the 
present  work. 

Failure  of  Rocks  Under  Dynamic  Tensile  Loads 

High  rate  loading  occurs,  for  example,  when  explosive  in  contact 
with  a  sample  is  detonated  or  when  a  sample  is  subjected  to  high  velocity 
impact.  Under  these  conditions  fracture  is  induced  by  tensile  stresses 
arising  from  intersections  of  reflected  shock  waves.  Such  tensile  stresses 
typically  have  rise  times  of  the  order  of  tens  of  nanoseconds  and  endure 
for  several  tenths  of  a  microsecond,  during  which  time  activated  cracks 

grow.  The  extent  to  which  individual  cracks  propagate  and  coalesce  depends 

t: 

predominantly  on  the  duration  of  the  stress  pulse;  thus  time  as  well  as 
stress  plays  an  important  role  in  the  dynamic  fracture  criterion. 


3 


Dynamic  tensile  failure  in  metals,  where  large  numbers  of  micro¬ 
cracks  form  but  do  not  fragment  the  specimen,  has  received  considerable 

attention,  and  several  schemes  have  been  proposed  to  include  time  in 

7—13 

dynamic  fracture  criteria.  A  general  method  was  presented  by  Tuler 

and  Butcher14  and  by  Gilman  and  Tuler,15  proposing  that  a  damage  function 
<|>is  a  function  of  the  entire  stress  history 

.  t 

$=  J*  f(cr(t))  dt  (l) 

—00 

where  <J>  is  selected  to  be  any  convenient  function  of  damage.  For  example, 
0  could  be  the  total  number  or  volume  of  the  microcracks  formed,  It  was 
suggested  that  f(cr(tj)  be  expanded  in  powers  of  <J-Qq>  where  is  a 
threshold  stress  below  which  no  damage  will  occur  regardless  of  stress 
duration.  It  was  further  suggested  as  a  zero  order  approximation  that 
one  term  might  be  expected  to  be  dominant,  so  that 

t  x 

<!>  =  ./  <cr-xro)  dt 


where  \  and  a  are  considered  material  constants  determined  by  experiment, 
o 

Much  effort  has  been  expended  in  trying  to  relate  observed  shock-induced 

damage  thresholds  in  metals  to  Eq.  (2).  Usually  4>  is  picked  to  be  a  level 

of  microscopic  damage  visible  at  some  level  of  magnification,  and  \  is 

10,16-18 

chosen  to  be  a  number  that  fits  the  data.  The  results  have  been 

only  partially  successful,  which  is  perhaps  to  be  expected  since  the 
definition  of  $  is  very  qualitative  and  a  given  material  is  described 
by  only  two  parameters.  In  short,  Eq.  (1)  is  undoubtedly  general  enough 
to  be  correct  in  principle,  but  Eq,  (2)  does  not  contain  enough  information 
about  a  given  material  to  be  useful  other  than  as  a  method  for  describing 
experimental  results. 
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Recently  an  approach  was  introduced  at  SRI  that  is  consistent 

with  the  basic  approach  of  Tuler  and  Butcher,  brt  concentrates  on  more 
precise  physical  definitions  of  the  <J)  and  f  functions  in  Eq.  (1)  than 
has  been  used  before.  This  approach  has  yielded  promising  results  for 
both  ductile  and  brittle  metals.  In  addition,  the  simple  physical  model 
underlying  the  approach  allows  dynamic  fracture  phenomena  to  be  related 
to  quasi-static  fracture  behavior,  and  thus  suggests  the  possibility  of 
eventually  understanding  the  nucleation  and  growth  of  microcracks  resulting 
in  fracture  at  all  rates  of  loading.  The  details  of  the  SRI  dynamic 
fracture  model  and  its  application  to  rocks  are  discussed  below. 

The  SRI  Dynamic  Fracture  Model 

Postfracture  examination  of  metallic  and  polymeric  specimens 
subjected  to  flat  plate  impact  revealed  that  under  dynamic  loading 

19—22 

conditions  fracture  originates  at  many  sites  throughout  the  specimen. 

Each  small  crack  then  propagates  and  achieves  a  certain  size  before  the 
stress  pulse  vanishes.  We  noted  the  analogy  to  a  metallurgical  phase 
transformation,  in  which  second-phase  particles  nucleate  and  grow  within 
a  matrix  material  under  the  influence  of  a  thermodynamic  driving  force, 
and  developed  a  model  based  on  this  analogy.  Thus  the  development  of 
fracture  damage  was  treated  as  a  process  of  crack  nucleation  and  growth. 

We  have  not  yet  treated  in  detail  the  third  stage  in  the  fracture  mechanism, 
coalescence,  in  which  individual  cracks  meet  and  join  up.  This  3tage  is 
particularly  important  for  more  brittle  materials  in  which  fragmentation 
occurs. 

Stress  dependent  crack  nucleation  and  growth  rate  functions  have  been 

19—21  22 

experimentally  determined  for  several  metals  and  for  one  polymer. 

These  functions  govern  the  development  of  fracture  damage  and  can  be  used 
to  predict  the  extent  and  location  of  damage  that  will  result  from  an 
arbitrary  loading  ni story. 
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The  nucleation  rate  and  growth  rate  functions  were  determined  from 

experiments  using  known  stress  histories.  The  stress  history  was 

23 

calculated  using  a  PUFF  code  with  appropriate  constitutive  relations. 
Unmodified  constitutive  relations  for  the  material  can  be  used  to 
calculate  the  stress  history  for  small  amounts  of  damage.  However,  for 
even  moderate  damage  levels  the  effect  of  the  damage  on  the  stress  history 
is  significant,  and  a  modified  equation  of  state  should  be  used  that  treats 
the  material  as  a  porous  solid  consisting  of  damaged  regions  and  undamaged 
regions. 

After  loading  by  a  known  pulse,  the  specimen  is  sectioned  on  a 
plane  parallel  to  the  impact  direction  to  reveal  the  damage  in  cross  section. 
Ductile  fracture  damage  appears  as  spherical  voids,  and  brittle  fracture 
damage  appears  as  plane  cracks.  An  important  aspect  of  the  analysis  is 
determination  of  the  volume  distribution  of  damage  from  such  surface 
observations.  The  procedure  for  developing  a  predictive  capability  for 
fracture  by  stress  waves  was: 

•  load  dynamically  and  recover  specimens 

•  achieve  experimental  control  so  that  the  damage 
can  be  stopped  in  different  stages  of  growth 

•  describe  quantitatively  the  size  and  the  spatial 
distribution  of  damage  in  the  volume  of  the  material 

•  specify  the  macroscopic  stress  and  the  stress  duration 
at  any  location  in  the  specimen  under  conditions  of 
nucleation  and  growth  of  damage. 

Nucleation  rate  and  growth  rate  functions  have  been  determined 

completely  by  combining  the  information  gained  in  the  steps  above. 

Moreover,  theoretical  work  in  the  nucleation  and  growth  processes  has 

made  it  possible  to  obtain  estimates  of  these  functions  for  materials 

other  than  those  tested . 
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For  material  that  fractures  in  a  ductile  manner  under  dynamic 

loading,  the  experimental  data  suggested  that  the  void  nucleation 

•  20 
rate  N  has  the  following  form 


exp 


wrere  or  is  the  applied  dynamic  stress,  or  is  the  threshold  stress  for 

No 

void  nucleatior,  and  N  and  or,  are  constants.  This  form  of  nucleation 

function  is  consistent  with  static  results  in  which  no  voids  appear 

for  stresses  less  than  the  yield  or  ultimate  strength.  It  not  only 

effectively  described  the  nucleation  of  voids  in  ductile  materials, 

but  also  was  successful  in  describing  crack  nucleation  in  more  brittle 
20 

materials. 

The  growth  rsite  of  brittle  cracks  might  be  expected  to  differ 
considerably  from  that  of  a  void  in  a  ductile  material.  For  ductile 
metals  ,  voids  were  found  experimentally  to  expand  according  to  a  viscous 
growth  law 


where  <t  is  the  threshold  stress  for  void  growth  and  T]  is  the 

Go 

material  viscosity.  Surprisingly,  materials  that  failed  by  propagation 
of  brittle-type  cracks  were  also  well  described  by  this  relation  and 
were  poorly  described  by  brittle  growth  laws  based  on  assumptions  of 
elasticity. 

Application  of  the  SRI  Dynamic  Fracture  Approach  to  Rocks 

In  any  attempt  to  model  dynamic  fracture  of  rock  materials,  several 
new  features  must  be  taken  into  account.  First,  a  much  more  pronounced 
inherent  defect  structure  exists  in  rocks  than  in  fabricated  materials 
such  as  metals  and  polymers.  The  defect  structure  profoundly  influences 
the  strength  und  fragmentation  characteristics  of  the  rock  and  therefore 
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must  be  determined.  Second,  rocks  generally  fracture  in  a  more  brittle 
manner  with  more  rapid  crack  velocities.  Thus  it  may  be  difficult  to 
stop  individual  cracks  from  growing  together  and  extending  to  specimen 
boundaries,  making  it  difficult  to  recover  intact  specimens  for  damage 
evaluations.  Finally,  increased  propensity  of  cracks  in  rock  to  run 
together  requires  not  only  that  crack  nucleation  and  growth  functions  be 
obtained,  but  also  that  the  third  stage  of  the  fracture  process,  crack 
coalescence,  be  treated.  A  correct  model  of  the  coalescence  behavior 
would  permit  calculation  of  the  fragment  size  distributions  resulting 
from  known  dynamic  tensile  loads. 

We  have  begun  in  the  present  work  to  apply  the  concepts  of  crack 
nucleation  and  growth  to  treat  dynamic  tensile  failure  of  rock..  The  steps 
taken  during  the  first  year  of  the  three-year  program  to  develop  a 
computational  model  for  high  rate  fracture  in  rock  were  as  follows: 

•  A  specimen  material  was  chosen  and  quantitatively 
characterized  with  respect  to  structure  and  properties. 

•  High  rate  tensile  tests  were  performed  to  measure  dynamic 
stress  histories  and  to  observe  the  fracture  behavior. 

•  Fractographic  analyses  were  carried  out  on  tested 
specimens  to  determine  the  fracture  mechanism. 

•  Construction  of  a  computational  model  was  begun 
for  predicting  the  failure  behavior  under  high  rate 
tensile  loads. 
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II  CHARACTERIZATION  OF  ARKANSAS  NOVACULITE 


ftlicrostructure 

Arkansas  novaculite— a  naturally  occurring,  polycrystalline  quartz-- 
is  a  mineralogically  simple  rock,  consisting  of  a  chemically  simple 
mineral,  and  was  therefore  a  suitable  material  with  which  to  begin  these 
studies.  Metallographic  examination  of  polished  and  etched  surfaces  at 
1000X  shows  that  novaculite  consists  of  equisized,  equiaxed,  and  randomly 
oriented  quartz  grains  having  an  average  diameter  of  about  10  p., 

Figure  1  shows  a  polished  surface  that  was  etched  for  2  minutes  in 
40%  HF  at  room  temperature  to  reveal  the  grain  structure.  Black  areas 
are  holes  where  natural  flaws  intersect  the  surface  or  where  grains  have 
been  pulled  out  during  the  polishing  process.  Bright  areas  are  caused  by 
reflected  light  from  internal  surfaces.  Observation  of  fracture  surfaces 
with  the  scanning  electron  microscope  reveals  the  grain  structure  even 
more  clearly  (Section  IV).  Crack  propagation  occurs  intergranularly , 
exposing  individual  grains  in  the  fracture  surfaces, 

Defect  Structure 

From  casual  observation  of  polished  surfaces  with  the  vinaided  eye, 
it  is  apparent  that  natural  novaculite  has  a  defect  structure  and  that 
the  defects  are  strongly  oriented.  Figure  2,  a  composite  photograph  at 
low  magnification,  illustrates  this.  Closer  examination  with  a  microscope 
shows  that  defects  are  present  on  essentially  two  distinguishable  size 
levels:  the  larger  defects  have  dimensions  of  a  few  hundred  microns; 
smaller  defects  are  roughly  the  size  of  the  grains. 

By  exploiting  the  translucency  of  novaculitt  and  focusing  into  the 
material  to  a  depth  of  about  100  u,  we  found  thal  the  larger  flaws  exist 
in  two  dominant  shapes:  plates  and  rods.  Best  results  were  obtained  by 
viewing  the  specimen  in  reflected  polarized  light  through  a  microscope 
slide  and  a  film  of  matching  refractive  index  oil  (n  =  1.55)  on  the 
specimen  surface.  Figure  3  shows  a  rodlike  and  a  platelike  flaw  slightly 


9 


FIGURE  1  MICROSTRUCTURE  OF  ARKANSAS  NOVACUUTE  SHOWING  THE  EQUIAXED 
QUARTZ  GRAINS 
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500  n 

MP-1087-12A 


FIGURE  2  COMPOSITE  MICROGRAPH  OF  A  BLOCK  OF  NOVACULITE 
SHOWING  THE  PREFERRED  ORIENTATION  OF  THE 
INHERENT  DEFECTS 
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below  the  surface  of  polish  and  inclined  at  an  angle  to  it,  so  that  only 

a  section  of  each  is  in  focus.  The  planes  of  the  rather  homogeneously 

distributed  platelike  flaws  are  roughly  parallel  to  one  another,  and 

most  of  the  rodlike  defects  are  Inclined  at  about  45°  to  these  planes. 

(Rodlike  flaws  possibly  originate  from  tiny  bubbles  of  water  that  are 

rejected  during  crystallization  of  the  rock  and  coalesce  to  form  capillary 
24 

tubes.)  The  size  distribution  of  these  flaws  was  determined  and  is 
discussed  later.  The  plates  would  appear  to  be  particularly  effective 
crack  starters,  and  their  strong  preferred  orientation  ic  expected  to 
result  in  highly  anisotropic  fracture  properties. 

Determination  of  Maximum  Crack  Size 

As  has  been  mentioned,  the  size  of  the  largest  flaw  existing  in  a 
specimen  determines  to  a  large  extent  its  quasi-static  strength.  Thus 
in  establishing  a  failure  criterion  it  is  necessary  to  be  able  to 
establish  the  maximum  crack  size.  For  opaque  materials  this  is  a  formidable 
task.  A  number  of  techniques  have  potential  for  providing  pictures  of 
Internal  artifacts  in  opaque  materials.  These  include  methods  based  on 
holography,  ultrasonic  imaging,  infrared  photography  and  x-ray  topography. 
However,  these  imaging  methods  are  still  in  the  early  stages  of  development, 
and  it  is  difficult  to  have  much  confidence  in  the  accuracy  of  the  results. 

A  better  established  approach  is  to  section  the  rock  specimen  perhaps 
on  three  perpendicular  surfaces  and  measure  the  crack  traces  intersecting 
a  sufficiently  large  area  of  the  surfaces  of  section.  If  many  crack  traces 
are  observed  whose  average  length  is  not  greatly  exceeded  by  any  one  trace, 
then  statistically,  the  maximum  crack  diameter  is  approximately  equal  to 
the  length  of  the  largest  trace. 

Such  condi tiona  seem  to  prevail  in  Arkansas  novaculite;  seldom  was 
any  crack  trace  obaerved  whose  length  oxceeded  greatly  those  of  many  other 
crack  traces.  The  largest  penny-shaped  crack  in  the  sample  of  novaculite 
used  in  this  work  waa  established  in  this  way  to  be  about  500  ^  in  diameter. 
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Determination  of  the  Initial  Crack  Size  Distribution 


Under  high  rate  loading  conditions,  cracks  smaller  than  the  largest 
crack  may. be  activated  and  contribute  to  specimen  failure.  Thus  in 
developing  a  capability  to  predict  fracture  behavior,  of  rock  under 
dynamic  loads,  it  is  necessary  to  be  able  to  detexmine  the  size  distri¬ 
bution  of  Inherent  flaws.  As  was  done  in  the  determination  of  the  maximum 

crack  size,  novel  imaging  techniques  were  rejected  in  favor  of  the  cheaper 

and  ^ore  reliable  technique  of  statistically  transforming  artifacts  on 
surfaces  of  section. 

The  traces  of  cracks  as  seen  in  a  composite  photograph  of  a  sectioned 

and  polished  rock  surface  were  counted  and  measured.  These  data  were  then 

statistically  transformed  to  obtain  the  size  distribution  of  Inherent 
cracks  per  unit  volume. 

Nino  overlapping  photographs  at  40X  were  required  to  span  the  diameter 

of  a  rock  specimen.  The  total  photographed  area  was  slightly  more  than 
2 

1.1  cm  ,  and  in  this  area  194  preexisting  crack  traces  were  counted  and 

measured.  These  data  were  converted  by  means  of  a  Schell  type  statistical 
3,4 

transformation  implemented  by  the  BABS2  computer  code  developed  at  SRI; 

the  results  are  presented  in  Figure  4. 

Here  the  cumulative  concentration  of  cracks  having  radii  greater  than 

radius  is  plotted  as  a  function  of  c^.  The  relatively  few  (194)  traces 

2 

observed  on  about  1  cm  of  the  surface  of  section  transform  into  a  very 

3 

large  volume  density  (~100, 000/cm  ).  The  size  distribution  of  preexisting 
cracks  has  a  parabolic  form  in  log-normal  space  with  a  cutoff  at  about 
c  s  500  M.  The  curve  is  well  described  by  the  analytical  expression 

N  =  6(c)  expCll.l  -  (3.7  x  102)c  +  (0.42  x  104)c2]  (3) 

where  6(c)  a  1  for  0  <  c  £  0.05 
6(c)  «  0  for  c  >  0.05 
and  c  is  in  centimeters. 


14 


CUMULATIVE  NUMBER  OF  CRACKS  OF  RADIUS  GREATER  THAN  c  —  numtoer/cm 


I 


i 


Quasi -Static  Strength  Measurements 


*  ,  \ 

i  The  tensile  strength  of  novaculite  under  quasi-static  loading 

!  25,26 

conditions  was  determined  using  the  SRI  expanded  ring  test  shown 

1  t  s 

in  Figure  5.  In>  this  test  hydrostatic  pressure  acts  radially  against 

the  Internal  wall  of  a  cylindrical  specimen  to  create  a  uniform  tangential 

1 

tensile  stress  in  the 'specimen  wall.  Nonaxial  stresses  caused  by  mis- 

1  :  l 

alignment  and  localized  stress  concentrations,  which  normally  arise  from 
gripping  or,  supporting  the  test  specimens,  a*e  eliminated  in  this  method. 
Ring-shaped  specimens  were  diamond  ground  from  oversized  blanks  to 


the  following  dimensions : 


O'.D., 


Height 


2,000  ±  0.0005  in. 
2.300  *  0.0005  in. 

i 

0.300  ±  0.0005  in. 


-4  -1 

All  strength  measurements  were  made  at  a  strain  rate  of  3  x  10  sec 

1  '  l  i  -5 

Measurements  were  taken  both  in  air  and  under  a  vacuum  of  10  torr. 

'  '  27 

Some  materials  show  higher  strengths  under  vacuum,  indicating  that  stress 
corrosion  is  an  important  factor.  No : significant  difference  was  seen  for 

I  I  ,i 

novaculite.  Results  of  15  tests  in  air  indicated  a  tensile  strength  of 

2  2 
441  ±  20  kg/cm  and  12  tests  under  vacuum  indicated  448  ±  33  kg/cm  . 

1  2 

The  extremes  measured  were  517  and  372  kg/cm  .  These  strengths  are 

28 

considerably  higheij  than  those  reported  rby  Hardy  and  Jayaraman  for  a 

2  ' 

variety  of  rocks,  such  as  Barre  granite,  88  kg/cm  ;  Crab  Orchard  sandstone, 

2  2 
71  kg/cm  ;  and  Leuders  limestone,  47  kg/cm  , 


Inherent  Fracture  Resistance 

i  i  ' 

The  longitudinal  sound  speed  c  in  novaculite  was  determined  by  the 

i  i  5 

time-in-flight  ultrasonic  imethod  to  be  5.9  x  10  cm/sdc;  a  value  of 
3 

2.63)  g/cm  was  detejrmined  for  its  density.,  These  values  were  used  to 

5  2 

calculate  a  Young's  modulus  E.of  about  6.4  x  10  Kg/cm  .  Friedman 

1  i  29  1  1 

and  coworkers,  using  the  double  cantilever  beam  technique,  determined 

t 

the  fracture  surface  energy  in  Chilhowee  quartzite  to  be  about  5  ±  2  x  10 
ergs/cm2. 


i 


i 


i 


SCREEN 

SHIELD 


EXPANDABLE 
RUBBER  TUBE 


SPECIMEN 


FIGURE  5  THE  SRI  EXPANDED  RING  TEST 


The  intrinsic  fracture  resistance  of  the  rock  (yE)^,  which  relates 

8 

the  fracture  strength  to  the  maximum  crack  size,  is  then  about  2  x  10 
—3/2  -3/2 

dynes  cm  or  about  200  kg  cm  Substituting  this  value  into  the 

Griffith  equation 

a  =  [2yE/iru-i/)c  ]*  (4) 

crit  max 

where  v  is  Poisson's  ratio 


we  establish 


0  =  160  c 

crlt  max 


-1/2 


as  a  fracture  criterion  where  the  fracture  strength  a  is  in  kg/cm  and 

crit 

the  half  length  of  the  largest  crack  c  is  in  centimeters. 

max 

Another  way  to  evaluate  this  property  is  to  measure  the  size  of  the 
largest  crack  in  a  specimen,  determine  the  fracture  strength,  and  use 
fracture  mechanics  theory.  This  was  done  for  novaculite  for  the  quasi¬ 
static  case.  We  could  not  establish  fracture  origins  on  novaculite 
fracture  surfaces  and  were  therefore  unable  to  determine  if  edge  flt»vs  or 
internal  flaws  controlled  the  strength.  Thus  we  could  only  estimate  the 

fracture  toughness  K  of  novaculite  from  the  measured  fracture  strength 

Ic 

(t  and  the  half  length  of  the  largest  plate-like  flaw  c  ,  by  using 
f  max 

Sneddon's  relation  for  an  internal  penny-shaped  crack  in  an  infinite 

30 

medium  subjected  to  uniform  tension 


2  3/2 

K  =  -T-  <T 9!  110  kg/cm  (5) 

Ic  Vff  f  max 


as  a  lower  bound  and  the  expression  for  an  edge  crack  in  a  semiinfinite 

31 

sheet  subjected  to  tension 


K  *  1.12  a  Vffc  =*  175  kg/cm 
Ic  f  max 


3/2 


(6) 


as  an  upper  bound. 

32 

These  values  are  considerably  lower  than  the  value  Bieniawski 

3/2 

obtained  from  4-point  bend  tests  on  quartzite,  630  kg/cm  ,  or  that 

29 

calculated  from  Friedman  et  al  measurements  of  the  fracture  surface 

3/2 

energy,  280  ±  110  kg/cm  .  In  view  of  the  results  of  the  dynamic  tensile 

strength  measurements  reported  in  the  next  section,  we  feel  that  Eq.  (3) 

is  more  appropriate  than  Eq.  (4)  for  these  tests  and  hence  that  the  value 
3/2 

of  110  kg/ cm  is  the  more  accurate.  Thus  solving  Eq.  (3)  for  the  fracture 
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stress  a  we  obtain 
f 


3/2  -1/2 

Of  -  198  kg/cm  > 

for  the  quasi-static  fracture  criterion  for  novaculite. 

The  results  of  the  characterization  work  for  Arkansas  novaculite 

are  summarized  in  Table  I.  It  remains  now  to  determine  the  rate 

sensitivity.  For  rocks  having  low  rate  sensitivity,  quasi-static  fracture 

theory  may  describe  dynamic  behavior  sufficiently  well.  Rate  sensitive 

rocks  will  require  a  modified  relation.  A  relation  based  on  fracture 

mechanics  concepts  and  modified  to  include  strain  rate  effects  has  been 

33 

proposed  by  Birkimer. 


Table  I 


STRUCTURE  AND  PROPERTIES  OF  ARKANSAS  NOVACULITE 


Microstructure 

Material 

Grains 

Defect  structure 
Flaw  shape 
Size  range 
Size  distribution 
Maximum  plate  size 
Preferred  orientation 
Physical  properties 

Longitudinal  sound  speed 
Density 

Elastic  modulus 
Mechanical  properties 


High  purity  polycrystalline  quartz 
Equisized,  equiaxed,  about  10  |X  in  diameter 

Channels  and  plates;  grain  boundary  seams 
From  ~10  to  ^600 
Exponential  (see  Figure  4) 

~  500 

Very  high;  plates  oriented  on  parallel  planes 

5 

5.9  x  10  cm/ sec 
3 

2.63  g/Cm 

5  2 

6.4  x  10  kg/cm 


Quasi-static  fracture  strength  440  ±  20  kg/cm 
a 


Surface  energy 


4  ,  2 

5  ±  2  x  10  ergs/cm 


Inherent  fracture  resistance  ~110  kg/cm 


3/2 


Reference  29. 
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Ill  DYNAMIC  EXPERIMENTS 


a  ■ 


High  Rate  Loading 

High  rate  loading  of  rock  specimens  was  accomplished  in  a  gas  gun. 
Projectiles,  12  Inches  long  by  2\  Inches  in  diameter  for  the  smaller  gas 
gun  and  30  Inches  long  by  4  inches  in  diameter  for  the  larger,  were 
accelerated  down  the  evacuated  barrel  of  the  gun  by  the  sudden  release 
of  pressure  from  an  adjacent  pressurized  chamber  of  helium.  Flat  target 
specimens  were  impacted  with  thin  polyethylene  or  plexiglass  flyer  plates 
attached  to  the  leading  edge  of  projectiles  (Figure  6).  Under  such  flat 
plate  impact,  compressive  waves  initially  run  into  the  specimen  and 
projectile  head  to  produce  a  state  of  one-dimensional  compressive  strain 
Tension  is  produced  in  the  specimens  when  release  waves  running  inward 
from  the  free  surface  of  the  specimen,  or  from  the  boundary  of  the  lower 
shock  impedance  backing  material,  meet  similar  release  waves  running 
inward  from  the  back  surface  of  the  flyer  plate  (Figure  7).  Knowledge 
of  the  shock  impedance  (product  of  density  and  shock  velocity)  of  the 
rock  and  flyer  plate  (and  backing  material  when  used)  allows  the  flyer 
thickness  to  be  designed  so  that  tension  occurs  first  at  the  midplane  of 
the  specimen.  A  sufficiently  large  diameter-to-thickness  ratio  for  the 
specimen  ensures  that  the  tensile  strain  is  one-dimensional  by  preventing  un¬ 
loading  waves  from  the  specimen  periphery  from  reaching  the  interior 
during  the  tension  phase. 

The  target  specimens  were  short  cylinders,  usually  about  1/4-inch 
thick  and  1/2  to  1-1/2-inches  in  diameter.  Their  axes  were  carefully 
aligned  to  coincide  with  that  of  the  gun  barrel  to  ensure  flat  plate 
impact.  The  specimens  were  fit  tightly  into  constraining  rings  of 
aluminum  which  is  a  very  close  match  to  novaculite  in  shock  impedance. 

Most  of  the  aluminum  constraining  rings  had  an  8°  taper  on  the  outer 
circumference  and  were  press-fitted  into  a  larger  aluminum  plate 
serving  as  a  specimen  holder.  Others  had  no  taper  and  were  lightly  held 
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FIGURE  6  EXPERIMENTAL  ARRANGEMENT  FOR  HIGH  RATE  TENSILE  TESTING 
OF  CYLINDRICAL  ROCK  SPECIMENS  IN  A  GAS  GUN 
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FIGURE  7  SCHEMATIC  DEPICTION  OF  THE  STRESS  HISTORY  OP  A  SPECIMEN, 
FLYER  PLATE,  AND  BACKING  MATERIAL  CAUSED  BY  IK'PACT 
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in  the  specimen  holder  by  several  dab*  of  epoxy.  Upon  impact  with  the 
flyer  plate,  the  specimen  and  constraining  ring  fly  free  froai  the  aluminum 
holder  and  into  a  catcher  tank,  which  is  lined  with  rags  to  prevent 
subsequent  impacts  and  possible  uncontrolled  dsasge.  The  projectile  snd 
projectile  head  sre  prevented  from  entering  the  catcher  tank  by  the  steel 
piste. 

Stress  History  Measurement 

34,35 

Plezoresistsnt  ytterbium  gages,  recently  developed  at  Sill  to 

measure  compressive  stress  in  the  sub-  snd  low-kllobsr  rsnge,  were  used 
to  measure  the  stress  histories  (peak  v.tress  snd  stress  duration)  pro¬ 
duced  in  the  specimens  by  the  impact.  A  small  grid-like  gage  having  a 
resistance  of  1  ohm,  a  thickness  of  0.002  inch,  and  covering  a  0.2-  uy 
0.2-inch  area  (Figure  8)  was  cemented  to  the  back  surfaces  of  rock 
specimens  and  backed  by  a  plexiglass  plate.  Because  of  the  location  of 
the  stress  gages,  the  stresses  measured  are  those  in  the  plexiglass  and 
sre  compressive.  The  tensile  stress  in  the  rock  was  then  calculated 
from  the  relative  shock  impedances  of  specimen  and  backing  material.  An 
indirect  method  of  tenalle  atreaa  measurement  is  necessary  since  in- 
materlsl  gages  that  measure  tenile  stress  directly  do  not  exist. 

The  gage  record  for  a  specimen  in  which  spsll  fracture  has  occurred 
yields  valuable  information  in  addition  to  the  peak  stress  and  stress 
duration.  Recompression  waves,  emitted  as  spsll  crscks  form,  Impinge 
on  and  reload  the  gage  to  an  extent  proportional  to  the  dynamic  tensile 
strength  of  the  rock.  Furthermore,  the  slope  of  the  reloading  pulse  is 
an  indication  of  the  rate  of  fracture— a  sharp  rise  corresponding  to 
brittle  behavior.  Finally,  the  measured  stress  history  provides  Important 
information  about  the  constitutive  relations  for  the  rock  and  their  rate 
dependence . 
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Dynamic  Tensile  Strength  Measurements 

We  performed  43  tensile  experiments  on  Arkansas  novacullte  st  nigh 
loading  rates,  using  the  gss  guns.  Fourteen  of  these  were  instrumented 
with  ytterbium  stress  gages  to  determine  stress  histories  and  to  measure 
the  magnitudes  of  recompression  waves  produced  by  fracture.  An  additional 
20  were  uninstrumented  experiments  whose  goals  were  to  determine  the 
dynamic  fracture  strength  of  novacullte  and  sny  effects  of  defect  orientation 
on  fracture  strength.  The  remaining  nine  experiments  were  either  of  a  pre¬ 
liminary  nature,  to  determine  in  a  rough  way  the  Impact  velocities  and  other 
experimental  conditions  required  to  carry  out  the  proposed  program,  or  of 
an  investigative  nature  to  determine,  for  example,  the  reason  for  the 
occurrence  of  undesirable  rsdial  cracking. 

Instrumented  Experiments 

A  summary  of  the  14  experiments  instrumented  with  ytterbium  stress 
gages  is  given  in  Table  II.  The  first  5  experiments  were  carried  out  with 
the  2)-lnch  gas  gun;  the  rest  were  done  with  the  4-inch  gas  gun  to  accomodate 
larger  specimens  snd  thicker  annuli  in  an  attempt  to  eliminate  rsdial 
crscking.  In  esch  of  these  experiments  the  gage  was  used  as  s  back-surface 
gage  near  the  interface  between  the  novacullte  specimen  and  a  plexiglass 
plate. 

Triggering  snd  noise  problems,  respectively,  ruined  the  first  two 
attempts  st  dynamic  measurements  (Shots  4  and  5);  specimen  recovery, 
however,  wss  successful.  The  noise  problems , presumably  piezoelectric  in 
origin,  were  overcome  by  using  the  differential  recording  scheme  shown 
in  Figure  9  snd  by  sepsrstlng  the  gsge  from  the  rock  by  a  thin  (0.010 
inch)  plexiglass  lsyer.  With  these  improvements  s  good  record  wss 
obtained  for  shot  6  (Figure  10) . 
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FIGURE  10  EXPERIMENTAL  YTTERBIUM  GAGE  RECORD  FOR  SHQT  6 

!  ' 

A*  can  be  deduced  from  the  depicted  stress  history  In  Figure  7 

:  i 

the  gage  record  for  a  specimen  which  does  not  undergo  spall  fracture 

should  consist  of  three  peaks  followed  by  an  anomalous  stretching 

signal  as  the  gage  falls  In  tension.  Thus  the  peaks  recorded  for 

Shot  6,  Figure  10,  Indicate  that  no  spall  fracture  occurred,  a  conclusion  1 

which  was  verified  by  sectioning  the  specimen  and  finding  no  cracks.  1 

i 

Figure  11a  shows  the  gage  record  predicted  by  a  simple  computer  code  for  the 
case  when  no  fracture  occurs.  The  constitutive  relation  assumed  is  that 

i 

for  an  ideally  elastic  solid.  The  stress-time  profile  is  similar  to  the  i 
experimental  record  of  Figure  70;  the  slopes  and  rounding  of  'the  experi¬ 
mental  record  are  caused  by  tilt  (one  part  :of  the  gage  sees  the  peak  stress 

before  another  part),  and  by  inelastic  material  response.  According  to 

35  1 

our  present  calibration  the  peaks  on  the  record  correspond  to  80i,  100, 

2 

and  80  kg/cm  in  the  plexiglass  backing  material.  .  The  calculated  peaks 

2  i 

are  140,  90,  and  60  kg/cm  .  The  experimental  peaks  are  spaced  1.5  ^isec 

■  i 

apart, in  good  agreement  with  1.6  psec  for  the  computed  peaks. 

29  ,  i 


X  10° 


! 


Shot  11  at  about  the  same  stress  in  novaculite  gave  a  record  (Figure  12) 

similar  to  that  for  shot  6,  The  three  peaks  correspond  to  150,  220,  and  120 
,  2 

kg/cm  ,  respectively,  in  plexiglass.  However,  the  peaks  are  separated  by 
2.7  p,sec.  No  spall  cracks  were  found  in  this  specimen,  although  it  was 
damaged  by  radial  cracks  (see  Figure  21). 


FIGURE  12  YTTERBIUM  STRESS  GAGE  RECORD  FOR  SHOT  11 


Shot  12  was  impacted  at  higher  velocity  than  Shot  11  and  should  have 

,  2 

experienced  a  peak  tension  of  about  520  kg/cm  .  The  gage  record  shows  a 
jagged  peak  too  high  in  amplitude  and  too  broad  in  time.  Once  again,  no 
distinct  spall  cracks  were  found,  but  the  specimen  was  seriously  damaged 
by  radial  cracking.  (See  discussion  on  origin  of  radial  cracks  later  in 
this  section.) 

At  this  point  we  decided  that  we  would  have  to  reduce  the  radial 
cracking  to  make  meaningful  measurements  of  the  fracture  strength.  Several 
experiments  indicated  that  thicker  aluminum  retaining  rings  or  specimens 
having  a  larger  diameter-to-thickness  ratio  would  keep  the  novaculite 
specimen  in  one-dimensional  compression  long  enough  to  prevent  radial 
cracking  from  release  waves  originating  at  the  outer  ring  boundary.  The 
shock  impedance  of  aluminum  matches  that  of  the  novaculite  closely,  so 

we  did  not  expect  difficulties  at  the  interface. 
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The  next  series  of  experiments  with  gages  (Shots  15-23)  were  conducted 

on  the  4-inch  gas  gun  with  specimens  1.0  or  1.25  inch  in  diameter  and  0.3 

inch  thick.  These  had  aluminum  retaining  rings  with  0.6-inch  wall  thickness 
o 

md  8  taper.  Our  main  goal  in  these  experiments  was  to  record  the  spall 
signal.  Consequently  all  the  tests  were  at  projectile  velocities 
considerably  above  that  required  to  cause  spall  fracturing. 

When  spall  cracking  occurs,  recompression  waves  are  emitted  as  the 

■  j  n 

cracks  open  up.  In  this  case  the  gage  would  see  the  initial  compressive 
pulse  from  the  impact  followed  immediately  by  a  lower  compressive  pulse, 
the  spall  signal,  as  the  waves  formed  by  the  fracturing  material  reach 
the  gage.  Figure  lib  shows  the  approximate  wave  profile  predicted  by 
the  computer  for  the  case  where  spall  cracking  occurs. 

Seven  of  the  experimental  records  are  shown  in  Figures  13  through  19. 

The  best  of  these  are  considered  to  be  for  shots  15,  17,  22,  and  23.  There 

is  a  general  tendency  to  have  a  main  peak  followed  by  two  lesser  peaks,  which 

we  interpret  as  spall  signals.  For  Shots  15,  17,  and  22  the  lesser  peaks 

occur  at  about  1.6  and  2.5  (isec  after  the  main  peak. 

For  Shot  15  (Figure  13)  the  amplitude  of  the  first  peak  corresponds 
2 

to  530  kg/cm  in  plexiglass,  while  the  peak  stress  calculated  from  the 

2  ,  2 

impact  velocity  is  1070  kg/ cm  in  novaculite  and  370  kg/cm  in  plexiglass. 

The  amplitudes  of  the  two  following  peaks,  each  measured  froir.  its  fore- 

2 

going  trough,  correspond  to  380  and  440  kg/cm  in  novaculite,  in  good 
agreement  with  the  measured  dynamic  tensile  strength  of  novaculite 
reported  in  the  next  subsection. 

The  record  for  Shot  16  (Figure  14)  is  noisy  and  the  signal  arrived 
13  g,sec  later  than  expected.  A  smaller  peak  does,  however,  follow  the 
main  peak  and  may  be  a  spall  signal. 
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The  record  for  Shot  17  (Figure  15)  is  similar  to  that  for  Shot  15 


The  main  peak  appears  too  high  (560  kg/cm  in  plexiglass,  while  calculation 

2 

from  impact  velocity  gives  320  kg/cm  in  plexiglass),  but  the  two  following 

2 

peaks  appear  reasonable  and  indicate  290  and  200  kg/cm  in  novaculite, 
measuring  from  the  foregoing  trough  to  the  peak  in  each  case. 


FIGURE  15  YTTERBIUM  STRESS  GAGE  RECORD  FOR  SHOT  17 


The  gage  in  Shot  18  was  noisy,  and  we  made  no  attempt  to  interpret 

the  record.  In  Shot  19  no  record  was  obtained  because  the  oscilloscope  did 

not  trigger.  The  record  for  Shot  20  (Figure  16)  appears  smeared  out.  The 

2 

main  peak  corresponds  to  240  kg/cm  in  plexiglass,  the  peak  stress  in 

2 

plexiglass  calculated  from  impact  velocity  is  300  kg/cm  .  The  plateau 

beginning  at  1.1  psec  after  the  main  peak  and  the  high  peak  3  psec  after 

the  main  peak  may  be  spall  signals.  The  latter  would  correspond  to  660 
,  2 

kg/cm  in  novaculite. 
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FIGURE  16  YTTERBIUM  STRESS  GAGE  RECORD  FOR  SHOT  20 


The  record  for  Shot  21  (Figure  17)  Is  similar  to  that  of  Shot  20, 

2 

except  that  the  amplitude  of  the  main  peak  is  very  low,  70  kg/cm  in 

2 

plexiglass  as  compared  with  the  calculated  value  of  330  kg/cm  .  The 
same  sort  of  plateau,  sharp  trough,  and  major  rise  follow  the  main  peak 
as  in  Shot  20. 


FIGURE  17  YTTERBIUM  STRESS  GAGE  RECORD  FOR  SHOT  21 
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The  records  for  Shots  22  and  23  (Figures  18  and  19)  are  more  like  those 

for  Shots  15  and  17.  Two  smaller  peaks  follow  the  main  one.  The  timing  of  the 

peaks  for  Shot  22  agrees  with  Shots  15  and  17;  the  peaks  in  Shot  23  are  more 

2 

spread  out  in  time.  The  main  peak  for  Shot  22  corresponds  to  420  kg/cm 

2 

in  plexiglass  as  compared  with  330  kg/cm  calculated  from  impact  velocity. 

The  first  following  peak  is  nearly  flat,  and  the  second  corresponds  to 
2 

270  kg/cm  in  novaculite.  !?or  Shot  23,  the  main  peaks  agree  well;  the 

2 

gage  recorded  290  kg/cm  in  plexiglass  and  the  calculated  peak  is  270 

2  2 
kg/cm  .  The  first  following  peak  corresponds  to  150  kg/cm  and  the  second 

2 

to  420  kg/cm  in  novaculite. 

In  summary,  the  gage  records  clearly  show  the  difference  between 
the  stress  profile  when  no  spall  occurred,  such  as  Shot  6,  and  the  profile 
with  spall,  such  as  Shot  15.  A  number  of  difficulties  were  experienced 
in  obtaining  good  records.  Most  of  these,  we  feel,  were  technique  develop¬ 
ment  problems  that  have  been  overcome.  We  are  now  confident  that  the  gage 
is  recording  the  stress  profile  with  fair  accuracy.  However,  the  profile 
is  complicated  because  the  fracture  process  itself  is  complicated.  A 
number  of  spall  cracks  are  activated,  propagate,  and  coalesce. 

. Instead  of  a  single  spall  signal ,  which  would  be  expected  if  a  single 
spall  crack  were  formed,  we  see  several  spall  signals  superimposed  on 
each  other.  The  maximum  amplitude  of  a  spall  signal  on  good  records 
corresponds  to  the  fracture  strength,  as  expected.  In  future  experiments, 
with  our  recent  Improvements  in  the  specimen  release,  we  expect  more 
interpretable  results.  Shots  can  be  fired  with  peak  stresses  closer 
to  the  spall  threshold,  causing  less  fracture  damage,  and  consequently 
generating  a  simpler  stress  profile.  Furthermore,  the  computed  stress 
histories  for  the  above  shots  were  based  on  oversimplified  constitutive 
relations  that  did  not  contain  stress  relaxation  caused  by  the  developing 
damage.  Computations  with  more  sophisticated  constitutive  relations  .ire 
presently  underway. 
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Uninstrumented  Experiments 

A  series  of  20  uninstrumented  experiments  was  performed  to  measure 
the  dynamic  tensile  strength  of  novaculite.  The  two  main  goals  were  to 
determine  the  effect  of  strain  rate  on  the  fracture  strength  of  novaculite 
and  to  determine  the  effect  of  defect  orientation  on  the  dynamic  fracture 
strength.  A  third  goal  was  to  produce  fracture  surfaces  for  examination 
in  the  scanning  electron  microscope  to  gain  Information  concerning  the 
fracture  mechanism. 

Ten  specimens  were  cut  so  that  the  Impact  direction  was  normal  to 

the  planes  of  the  initial  cracks;  crack  planes  in  the  remaining  ten  were 

oriented  so  as  to  contain  the  impact  direction.  All  specimens  were  1/2 

inch  in  diameter  by  1/4  inch  thick  and  press-fitted  into  1-1/2-lnch 

diameter  by  1/4-inch  aluminum  annuli .  The  outer  periphery  of  each 

o 

annulus  was  provided  with  an  8  taper  to  facilitate  ejection  upon  Impact 

of  the  specimen-ring  assembly  from  the  aluminum  target  plate. 

Very  extensive  cracking,  usually  resulting  in  blow-out  of  one  or 

both  sides  of  the  specimen  from  the  aluminum  annulus,  was  produced  in  sll 

2 

of  the  first  nine  experiments  even  at  stresses  as  low  as  200  kg/cm  . 

This  initially  puzzling  behavior  was  prevented  io  the  subsequent  11  shots 
by  lightly  tacking  the  tapered  annuli  to  the  targe*  plates  in  three  or 
four  places  with  epoxy,  as  explained  in  the  next  section,  Instead  of 
press-fitting.  In  the  latter  experiments  no  radial  cracking  occurred,  and 
a  good  estimate  of  the  dynamic  fracture  strength  was  obtained. 

The  results  of  this  series  of  experiments  are  presented  in  Table  III. 

The  dynaadc  tensile  strength  was  taken  as  the  average  of  the  highest  stress 
at  which  no  damage  could  be  observed  on  diametrical  sections  of  the 
specimen,  and  the  lowest  stress  at  which  incipient  spallation  occurred. 

A  cross  sectional  rlew  of  incipient  spallation  is  provided  in  Figure  20. 

2 

For  the  parallel  specimen  orientation,  specimen  36  did  not  crack  at  402  kg/cm  , 

2 

but  specimen  35  did  crack  at  444  kg/cm  ;  for  the  normal  specimen  orientation, 
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FIGURE  20  CROSS  SECTIONAL  VIEW  OF  SPECIMEN  39 
SHOWING  INCIPIENT  SPALLATION 

no  cracking  was  observed  In  specimen  43,  which  was  subjected  to  a  peak 

2 

tensile  stress  of  396  kg/cm  ,  whereas  specimen  37  showed  cracking  at 
2 

433  kg/cm  .  Thus  the  dynamic  tensile  strength  of  novacullte  In  the 

2 

direction  normal  to  the  planes  of  Inherent  cracks  Is  415  ±  20  kg/cra  ; 

2 

the  dynamic  strength  In  the  transverse  direction  Is  423  ±  20  kg/cm  . 

We  conclude  that  the  dynamic  tensile  strength  of  Arkansas  novacullte  Is 
Insensitive  to  the  orientation  of  preexisting  flaws.  In  view  of  th 
pronounced  crack  orientation  anisotropy  of  novacullte,  the  fracture 
strength  Isotropy  was  surprising.  The  probable  explanation,  however, 

Is  that  the  tips  of  preexisting  Intergranular  flaws  lie  at  high  angles 
to  the  resultant  macroscopic  fracture  plane  and  are  equally  difficult 
to  activate  In  the  two  orientations.  It  should  be  borne  lr,  mind  that  the 
stress  state  Is  one  of  uniaxial  strain,  and  thus  trlaxlal  tension.  This 
means  that  the  stresses  In  the  two  perpendicular  directions  normal  to  the 
Impact  direction  are  nearly  as  large  as  the  Impact  stress. 
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The  strain  r*ie  of  the  dynamic  tests  was  estimated  from  the  oscill¬ 
oscope  records  of  the  ytterbium  stress  gages  in  the  instrumented  experi¬ 
ments  described  above.  The  rise  time  of  the  stress  pulse  is  about  a 

-4 

tenth  of  a  microsecond,  the  elastic  strsin  is  about  7  x  10  ,  hence  the 

3  -1 

strsin  rate  is  estimated  to  have  been  approximately  7  x  10  sec  ,  since 
this  is  more  than  seven  orders  of  magnitude  higher  than  the  strain  rates 
of  the  quasi-static  ring  tests,  the  negligible  difference  in  strength 

values  indicates  that  the  fracture  strength  is  also  strain-rate  ln- 

a 

aensitive,  at  least  in  this  rsnge. 

Origin  of  Radial  Cracks 

In  the  early  stages  of  this  work  specimens  subjected  to  dynamic 
tensile  loading  in  the  gas  guns  were  observed  to  have  fractured  in  a 
radial  mode  in  addition  to  the  planned  spallation  mode.  Under  conditions 
of  flat  plate  impact,  the  cylindrical  rock  specimens  were  expected  to  fail 
by  propagation  of  so-called  apall  cracks  roughly  normal  to  the  impact 
direction.  It  was  therefore  somewhat  surprising  to  consistently  observe 
a  radial  network  of  cracks  parallel  to  the  Impact  direction  and  extending 
from  the  impact  side  into  about  the  midplane  of  the  specimen.  Figure  21 
is  a  photograph  of  the  radlsl  network  of  cracks  on  the  impact  surface  of 
a  specimen,  and  Figure  22  ia  a  diametrical  section  showing  the  depth  to 
which  the  cracks  extend  into  the  interior. 

Such  radial  cracks  were  undesirable  since  they  complicated  the  failure 
proceas  and  made  interpretation  of  stress  gage  records  uncertain.  The 
reason  for  their  occurrence  was  elusive.  We  first  speculated  that 
rarefaction  waves  running  in  from  the  specimen  periphery  were  responsible 


36 

In  light  of  experimental  resylta  Rinehart  on  limestone  and  marble 
at  strain  rates  exceeding  10  sec  ,  it  is  conceivable  that  novaculite 
could  show  rate  sensitivity  at  higher  strain  ratea. 
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MS-10S7-10 


FIGURE  21  RADIAL  CRACKING  PATTERN 
ON  THE  IMPACT  SURFACE 
OF  NOVACULITE 

for  radial  cracking,  so  in  attempting  to  eliminate  them,  specimens  were 
inserted  into  thick  aluminum  annuli.  Since  the  shock  impedance  of 
aluminum  is  similar  to  that  of  novaculite,  this  alteration  had  the 
effect  of  displacing  the  specimen  periphery  so  thot  the  tensile  pulse 
producing  the  spall  cracks  was  completed  before  the  rarefaction  waves 
could  reach  the  specimen  interior.  An  equivalent  measure  taken  to 
eliminate  radial  mode  cracking  was  to  perform  experiments  in  a  larger 
bore  gas  gun,  using  larger  specimens. 

Several  experiments  using  this  technique  were  successful  in  that 
apall  cracks  but  no  radial  cracks  were  produced.  In  other  cases,  however, 
the  converse  was  true,  and  it  was  obvious  that  the  problem  was  not  solved. 
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FIGURE  22  CROSS  SECTIONAL  VIEW 

OF  A  SPECIMEN  SHOWING  K 
RADIAL  CRACKS  EXTENDING 
INJO  ABOUT  THE  MlOPLANE 
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During  the  course  of  carrying  out  the  aeries  of  uninstrumented  recovery 
shots,  an  observation  was  made  that  indicated  the  solution.  An  impacted 
specimen  was  not  ejected  into  the  catcher  tank,  but  rather  remained  in 
the  target  plate.  Radial  cracking  was  significant  and  suggested  that  the 

press  fitting  prevented  the  specimen  from  immediately  popping  free  under 

,  1 

the  low, impact  velocities,  and  the  ejection  delay  permitted  bending  waves 
to  produce  the  raldial  cracks.  To  tbst  this  theory,  specimens  were  weakly 
held  in  the  target  plate  with  a  dab  of  epoxy  in  several  places.  This 
technique  was  used  for  the' last  11  specimens  with  complete  success;  i.e., 

^  i 

no  radial  'cracking  occurred  in  any  of  these  experiments  (see  for  example, 
Figure  20). 


i  i 


i 


i  i 
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IV  THE  MECHANISM  OF  FRACTURE 


Fractographlc  Observations 

The  fracture  damage  produced  during  most  of  the  Impact  experiments 
consisted  of  many  cracks  in  essentially  two  main  orientations:  cracks 
lying  perpendicular  to  the  impact  direction  (generally  referred  to  as 
spall  cracks*  and  unexpected  cracks  whose  planes  roughly  contained  the 
Impact  direction.  The  spurious  radial  cracks  that  occurred  because  the 
specimens  were  held  too  tightly  in  the  target  plate  were  eliminated  in 
later  experiments  as  explained  in  the  previous  section,  and  thus  an 
unambiguous  fractographlc  analysis  could  be  performed  on  the  spall  cracks. 

Fractography  of  specimens  that  remained  intact  after  Impact  was  per¬ 
formed  by  sectioning  the  specimens  on  a  diameter,  polishing  the  surface 
of  section,  and  examining  the  fracture  pattern  by  means  of  optical  micro¬ 
scopy.  At  low  stresses  Just  slightly  in  excess  of  the  dynamic  fracture 
strength,  only  a  few  spall  cracks  (usually  less  than  four  Intersected 
the  surface  of  section)  were  formed  (Figure  20) ,  but  these  were  usually 
large,  extending  typically  to  about  a  third  of  the  specimen  diameter. 

Some  evidence  of  crack  branching  was  observed,  indicating  that  propagation 
velocities  were  hlgh,probably  approaching  a  limiting  crack  velocity.  The 
low  number  of  cracks  was  expected  at  these  stresses,  since  only  the  largest 
or  near  largest  cracks  could  be  activated.  Furthermore,  the  stress 
quenching  effect  in  the  volume  around  the  rapidly  running  cracks  would 
tend  to  prevent  additional  near-critical  cracks  from  being  activated 
during  the  tensile  pulse. 

At  higher  stresses,  say  50%  in  excess  of  the  dynamic  tensile  strength, 
considerably  more  Inherent  cracks  were  activated.  Probably  growing  at 
near  limiting  velocities,  these  cracks,  because  of  their  high  number 
density,  soon  grew  into  each  other,  and  the  fracture  process  proceeded 
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to  the  next  stage— -coalescence.  Figure  23  is  a  photograph  of  a  polished 
cross  section  showing  coalesced  cracks  and  incipient  fragmentation. 

Individual  cracks  have  joined  up,  and  in  doing  so,  have  produced  a  net¬ 
work  of  weakly  held  unfractured  rock  fragments,  Ths  missing  half  of  the 
specimen  in  Figure  23  has  undergone  complete  fragmentation,  and  the  indi¬ 
vidual  fragments  have  fallen  out  of  the  restraining  ring. 

The  fracture  surfaces  of  spalled  and  fragmented  novaculite  specimens 
were  examined  by  optical  and  scanning  electron  microscopy.  The  purpose 
was  to  look  for  inhomogeneities,  which  could  have  served  as  crack  acti¬ 
vation  sites,  and  such  crack  growth  markings  as  river  lines  and  hesitation 
bands,  all  of  which  would  provide  clues  concerning  tho  fracture  mechanism. 

As  is  typical  of  most  rocks,  the  fracture  surfaces  vere  nearly 
featureless  and  yielded  scant  evidence  of  how  failure  occurred.  The 
broken  surfaces  consisted  of  countless,  well-defined,  equlsized  polygonal 
blocks — quartz  grains  exposed  by  t.he  passage  of  a  crack  along  the  grain 
boundaries  , as  illustrated  by  the  background  structure  in  Figure  24.  Thus 
the  measured  fracture  strength  is  really  the  strength  of  the  '  cement 
holding  individual  grains  together,  and  not  of  the  quartz  itself.  No 
obvious  defects  could  be  found  where  cracks  had  nucleated,  and  no  telltale 
lines  such  as  often  mark  the  boundaries  of  preexisting  flaws  in  other 
materials  when  activated  under  a  tensile  pulse  were  evident.  Although 
no  fracture  markings  could  be  found  to  confirm  it,  it  is  likely  that 
anomalously  large  grains  such  as  that  shown  in  Figure  24  served  as  acti¬ 
vation  sites. 

It  was  possible  on  occasion  to  identify  propagation  markings  and 

hesitation  lines  on  fracture  surfaces,  Figure  25,  which  suggest  that 

cracks  expand  radially  outward  from  their  activation  sites  on  planes 

normal  to  the  tensile  direction.  This  crack  propagation  mode  was  found 

22 

to  operate  in  polymers,  Figure  26,  under  similar  dynamic  loading  conditions. 
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FIGURE  24  ELECTRON  MICROGRAPH  OF  A  LARGE  QUARTZ 
GRAIN,  FOUND  ON  THE  FRACTURE  SURFACE 
OF  AN  IMPACTED  SPECIMEN,  THAT  MAY  HAVE 
SERVED  AS  A  SPALL  CRACK  NUCLEATION  SITE 
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FIGURE  25a  SECTION  OF  CIRCULAR  HESITATION  LINE  ON  THE  FRACTURE  SURFACE 
OF  DYNAMICALLY  LOADED  NOVACULITE 
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FIGURE  26b  SECTION  OF  CIRCULAR  HESITATION  LINE  ON  THE  FRACTURE  SURFACE 
OF  DYNAMICALLY  LOADED  NOVACULITE 
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FIGURE  28  INTERNAL  CRACKS  IN  POLYCARBONATE  PRODUCED 
BY  A  SHORT-LIVED  TENSILE  PULSE 
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Fracture  Mechanism  in  Novaculite  Under  a  Short  Tensile  Pulao 


On  the  basla  of  these  f ractographlc  observations  and  previous 
observationa  carried  out  on  a  transparent  polymer,  we  propose  that  frocturo 
in  rocks  under  high  rate  loading  occurs  in  the  four  stagoa  desci lbod 
below. 

1.  Simultaneous  activation  of  a  number  of  preexisting  structural 
defects — The m*  include  channel-like  pores,  large  voids,  planar  fissures, 
and  grain  boundary  aeama  (see  Section  III).  Since  the  stress  required 
to  activate  a  flaw  decreases  with  Increasing  flaw  size  and  decreasing 
root  radius  of  the  crack  periphery,  the  largest  and  sharpest  flaws  are 
expected  to  be  the  moat  influential  in  instigating  the  fracturo  process. 
However,  in  the  case  of  novaculite,  the  dynamic  strength  was  found  to  be 
Independent  of  the  orientation  of  the  large  plate-like  flaws,  Indicating 
that  the  randomly  orientated  grain  boundary  flaws  may  control  internal 
spall  initiation. 

2.  Growth  of  the  activated  flaws  radially  outward — Circular  hesi¬ 
tation  lines  on  fracture  surfaces  of  Impact  loaded  novaculite,  Figure  25, 

suggest  that  penny-shaped  cracks  grew  by  radial  expansion.  This  con- 

22 

elusion  ia  supported  by  similar  impact  experiments  on  a  transparent 
plastic  where  individual  cracks  could  be  directly  observed,  Figure  26. 

The  path  of  the  fracture  la  determined  by  two  factors:  the  direction  of 
the  maximum  tensile  stress  and  the  local  structural  discontinuities  in 
the  vicinity  of  the  crack  periphery.  Strain  energy  release  per  unit 
crack  advancement  is  greatest  for  crack  extension  normal  to  the  direction 
of  maximum  tension.  This  tends  to  maintain  crack  growth  on  a  level  plane 
since  the  stress  pulse  in  theae  experiments  la  constant  in  direction. 

The  absence  of  smooth  planar  cracks  in  novaculite  ia  a  result  of  the 
micro  and  defect  structure.  Since  the  strength  of  the  quartz  grains  is 
greater  than  the  strength  of  the  bonding  between  them,  the  microstructure 
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determines  the  small  acale  roughneaa  of  th*  fracture  path,  because  the 
advancing  crack  must  paaa  around  the  grains  rather  than  through  then. 

In  a  alnllar  manner  the  defect  structure  cauaea  both  small  and  large 
scale  deviation  fjrom  the  planar  fracture  path  favored  by  the  maximum 
tenallo  atreaa.  The  crack  intersects  inherent  pores,  fissures,  and  voids 
aa  It  grows,  and  incorporates  these  into  the  crack.  Such  flaws  do  not 
generally  llo  In  the  tensile  plane,  however,  and  thus  roault  In  larger 
scale  Irregularity  of  the  fracture  surface.  Moreover,  since  the  crack 
prefers  to  follow  a  leaat-roslstanco  path,  it  Is  not  even  necessary  that 
structural  discontinuities  Intersect  the  plane  of  preferred  propagation; 
rather  It  la  often  sufficient  when  a  large  flaw  lies  somewhat  out  of  the 
preferred  plane,  for  an  interaction  of  the  stress  fields  can  cause  the 
crack  to  turn  into  the  flaw.  Considerable  roughness  also  results  from  a 
apllttlng-up  of  the  main  crack  Into  two  or  more  branches  as  a  consequence 
of  the  defect  structure  or  of  high  fracture  spoeds.  Secondary  branches 
may  run  a  distance  and  atop,  or  they  may  Join  up  with  other  cracks,  which 
sometimes  produces  free  chips  of  material.  Crack  branching  Is  responsible 
for  the  shattering  phenomenon  in  brittle  materials  under  quasi  static 
loading. 

The  fracture  path  thon  la  determined  not  simply  by  the  direction  of 

tho  maximum  tensile  stress  o  but  rather  by  the  ratio  of  o  to 

max  '  max 

the  local  strength  of  the  materlsl  S  ,  Fracture  propagation  occurs 

1  oc 

on  those  planes  that  are  normal  to  tho  direction  In  which  o  /s 

max  loc 

is  a  maximum.  Notice  that  since  varies  strongly  with  position  In 

the  material,  tho  fracture  path  will  not  be  constant  but  will  wander 
haphazardly  through  the  specimen  along  the  course  of  least  resistance. 

3.  Coalescence  of  the  expanding  cracks— Two  radially  expanding 
cracks  on  nearby  planes  will  eventually  approach  each  other  to  the  point 
where  the  concentrated  stress  fields  associated  with  the  crack  tips 
overlap.  The  material  between  the  two  cracks  is  subjected  to  magnified 
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itniNif  and  Imcium  of  tho  proximity  of  two  fro*  surfaces  (In  uhlch  the 
stress  haa  fallen  to  sero),  the  maximum  atroaa  direction  ia  altered. 
Accordingly,  the  direction  of  crack  extension  la  altered,  uaually  caualng 
one  or  both  cracka  to  turn  in  to  one  another.  Thla  proceaa,  carried  to 
completion  ao  that  the  two  cracka  join  up,  la  called  crack  coalescence 
and  la  an  Important  stage  in  the  fracture  process.  Inability  of  cracka 
to  coalesce  would  be  aianlfeated  by  a  significant  Increase  In  the 
toughness  of  rocks.  Coslescence  accelerates  fracture  since  ersek  slse 
dlscontlnuously  Increases,  resulting  in  a  sudden  reduction  In  the  avail* 
able  load-supporting  area  and  a  corresponding  step-llke  reduction  In  the 
material  strength.  The  crosa-sectlonal  area  capable  of  bearing  load 
diminishes  at  an  ever-accelerating  rate  and  leads  to  the  final  stage  of 
fracture. 

4.  Fragmentation  of  the  specimen— Tho  connectivity  of  Individual 
cracka  on  at  least  one  fracture  psth  la  coaiplete,  and  the  erseka  hsve 
run  to  the  outer  surfaces.  Crscka  on  other  planes  atop  growing  because 
the  stress  Is  relieved.  At  stresses  bsrely  exceeding  the  fracture 
strength,  the  specimen  stay  break  Into  only  two  pieces,  but  st  stresses 
considerably  greater  than  the  fracture  strength,  significant  fragmentation 
occurs.  The  musber  of  fragments  produced  Is  thought  to  be  largely  de¬ 
pendent  on  the  applied  stress  and  the  slse  distribution  of  lnltlsl 
defects.  Crack  branching  con  also  contribute  to  the  fragmentation  process. 
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V  A  COMPUTATIONAL  MODEL  FOR  DYNAMIC  FRACTURE  OF  ROCK 


For  cosquitatlonal  purposes  the  tensile  failure  of  rock  Is  assumed 
to  occur  in  four  stages: 

(1)  Activation  of  s  number  of  preexisting  structursl  defects 

(2)  Propagation  of  activated  cracks  radially  outward 

(3)  Coalescence  and  branching  of  propagating  cricks 

(4)  Isolation  of  Individual  rock  fragments  from  one  another. 

The  Important  quantities  to  be  determined  are) 

(1)  The  number  of  flaws  that  becoaie  running  cracks 

(2)  The  distance  that  each  crack  propagates 

(3)  The  degree  to  which  coalescence  and  branching  occurs 

(4)  Tbe  size  distribution  of  the  resulting  fragment  a. 

It  la  aasuwjd  that  theae  quantities  are  governed  by  three  factors:  the 
Initial  defect  structure  in  the  rock,  the  Inherent  fracture  resistance 
of  crack-free  material,  and  the  applied  stress  history.  To  predict 
the  failure  behavior  of  a  rock,  quantitative  aieaaurea  of  these  three 
factors  must  be  obtained  and  related  to  the  above  llstod  quantities, 

Which  are  then  related  to  each  other  by  meana  of  a  computational  model. 

During  this  year  the  first  two  stages  were  treated  In  some  detail; 
the  procedure  for  calculating  theae  quantities  and  the  results  for 
novacullte  are  presented  below.  The  latter  two  states  will  be  treated 
In  the  course  of  the  second  year. 

Stago  1:  Flaw  activation  calculations  for  novacullte — Vfe  assume 

flaw  activation  to  occur  according  to  Griff ith-Irwin  fracture  mechanics 
2-4 

concepts,  i.e„,  flaws  having  a  radius  c  larger  than  a  cortaln  size 

c*  *»  c* ( "0  are  activated  by  the  applied  stress  O  ,  whereas  smaller 

flaws  are  unaffected,  Frsctogrsphy  showed  that  internal  flaws  control 

the  activation  stage  of  the  fracture  process.  Of  the  various  geometries 

.  30 

for  which  fracture  mechanics  stresa  analyses  exist,  Sneddon  a  relation 
for  an  Internal  penny-shaped  crack  In  an  Infinite  elastic  medium  subjected 
to  a  uniform  tonslle  stress  normal  to  the  crack  plane 
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c*  * 


(7) 


/4 O2 
Ic 

■cst  nearly  applies  to  the  conditions  of  our  experiments,  l.e.,  the 

oracka  are  roughly  penny-shaped,  and  the  apeclnen  behaves  aa  if  it  were 

infinite  during  the  life  of  the  tenalle  pulse,  Using  this  expression  to 

relate  the  neaaured  dynamic  tenalle  strength  to  the  half  length  of  the 

largest  flaw,  we  oalculate  the  dynamic  fracture  toughness  of  novacullte 

3/2 

to  be  about  110  kg/cm  .  Substituting  this  value  into  Eq.  (7)  we  obtain 
the  functional  expression  for  relating  the  aize  of  the  smallest  pre  • 
existing  defect  which  will  be  activated  by  an  applied  dynamic  load  of 
stress  level  0 


c*  »  lO4/^  (8) 

To  obtain  the  number  of  defecta  which  will  be  activated  by  the 

applied  stress,  we  must  determine  the  aize  distribution  of  preexisting 

flaws  in  the  rock  aamplo  and  combine  this  information  with  Eq.  (7). 

Perhaps  the  most  reliable  and  convenient  wsy  to  determine  the  aize 

distribution  of  preexisting  flaws  la  by  sectioning  the  rock  sample, 

counting  and  measuring  the  flaw  tracea  on  the  section  surface,  and 

statistically  transforming  the  data.  Nine  overlapping  micrographs 

spanning  the  diameter  of  a  rock  specimen  were  assembled  into  a  composite 

picture  of  the  rock  surface.  The  total  area  was  slightly  greater  than 
2 

1.1  cm  and  in  this  area  194  tracea  of  preexisting  flaws  were  observed, 

counted,  and  measured.  We  assumed  that  crack  activation  occurs  most 

easily  at  the  large  plate-like  flawa,  although  in  light  of  the  constancy 

of  the  dynamic  fracture  strength  with  specimen  orientation  we  have  not 

ruled  out  tho  possibility  that  grain  boundary  fissures  may  control  the 

activation  stage  of  the  fracture  process.  The  data  were  converted  to 

obtain  the  size  distribution  per  unit  volume  by  means  of  a  statistical 

5,6 

transformation  similar  to  that  used  by  Schell.  The  procedure  was 
Implemented  by  a  computer  code. 
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The  cumulative  number  density  ol  flaws  haying  half  lengths  greater  , 
than  c  aa  a  function  of  c  are  presented  in  Figure  4  and  is  well 
described  by  the  analytical  form 

N  -  6(c)exp[ll-(3.7  x  102)c  +  (0,42  x  104)C2J  ! (4) 

'  ’  . 

where 

I 

6(c)  •-  1  for  0  <  c  ^  0,05  ■ 

6(c)  =  0  for  c  >  0,05 

•  i 

l 

and  c  la  In  centimeters. 

\  ! 

To  obtain  the  number  of  flaws  activated  N  by  an  arbitrary  i 

act 

dynamic  stress  o  ,  Eqs.  (4)  and  (8)  are  combined, 

I 

Nact  "  -  <3*7  x  +  (4.2  x  lO11)/^2]  (9) 

i  1  i 

l  , 

Here  , 

'  ♦ 

,  2  '  1 
6(0)  =  1  when  0  *  410  kg/cm 

i 

•  i 

and 

6(0)  *o  when  O  <  410  kg/cm2 

■  i 

Stage  2:  Crack  propagation  calculations  for  novacullte — The 

distance  each  crack  can  propagate  depends  on  the  crack  velocity  and  the 

♦  1  , 
duration  of  the  stress  pulse.  Ytterbium  stress  gages  were  used  to 

measure  the  latter  In  this  work  and  crack  velocities  were  Inferred  from  1 

direct  measurements  of  the  distance  propagated  by  cracks  In  essentially 

crack-free  material.  The  radii  of  the  cracks  shown  in  Figure  20  are 

i 

about  2  mm.  The  ytterbium  stress  gage  records  indicated  that  the  stress 

1  i 

duration  waa  about  one  microsecond,  In  agreement  with  the  result  obtained 

from  simple  calculations  using  the  measured  elastic  wave  speed  and  the 

thicknesses  of  the  specimen  and  flyer  plate  as  depicted  In  Figure  7. 

5 

Thua  crack  propagation  velocities  of  2  x  10  cm/sec  a (re  lndlcatedr-an 

( 

interesting  result,  since  It  Is  approximately  one-third  of  the  measured 

*  l 

l 

» 

:  ! 
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i 


longitudinal  wave  velocity  for  novacullte,  and  la  thus  in  agreement  with 

theoretical  estimates  of  the  maximum  crack  velocities  of  brittle  materi- 
,  37-40 

ais.  Materials  such  as  Armco  iron  and  beryllium,  investigated  in 

other  projects,15'41  fractured  in  a  brittle  manner  at  high  strain  rates, 
but  exhibited  viscous  crack  propagation  and  maximum  crack  speeds  well 
below  c^/3.  Knowing  the  crack  propagation, speed  and  that  propagation 

J 

occurs  radially  outward  from  activation  sites,  the  fracture  surface  area 

:  i 

produced  per  crack  A  during  the  period  of  the  stress  pulse  At  can 

i  1  i 

be  approximated  by 

i 

i  i- 


I  «  2fT(c{At/3)2  Q0) 

i  i 

i  1 

providing  the  crack  iwas  not,  stopped  prematurely  by  barriers  to  crack 

|  i 

growth  or  coalescence  with  other  cracks.  Thus  for  stresses  only 
,  slightly  greater  than  the  dynamic  tensile  strength  where  significant 

'  *  T 

crack  coalescence  do4s  not  occur,  the  total  fracture  ;area  produced 
in  novacullte  by  a  stress  pulse  of  knowii  peak  stress  a  and  duration 
At  is  given  by  . 

•  1  I  1 


i 


\ot:al 


2N 

act 

3 


(o^A,)2' 


(11) 


i  where  c^  is  the  longitudinal  wave  st>eed. 

i 

In  Eqq,  (10)  and  (11)  we  have  neglected  crack  acceleration,  assuming 

in  essence  that  cracks  reach  the  maximum  speed  c./3  in  zero  time.  A 
i  39  ^ 

solution  of  the  Dulaney  pnd  Brace  formula  confirms  that  the  acceleration 

.  *  •  . 

time  is  very  short  relative  to  the  stress  duration  in  this  work  and  thus 

*  I 

supports  this  assumption; 

t 

Thus  Eq.  (11)  uses  traditional  fracture  mechanics  and  simple  elastic 

.  •  -  ’  i 

crack  propagation  theory  to  obtain  a  crude  approximation  of  the  fracture 

i  .  ; 

*  . 

!  I  I 


1 


I 


damage  produced  In  novacullte  by  a  known  stress  history.  To  Illustrate 

«  2 

the  method  we  consider  a  square  tensile  pulse  500  kg/cm  in  amplitude  and 
0.8  microsecond  in  duration  and  calculate  step-by-step  the  amount  of 
fracture  surface  area  which  will  result.  According  to  Eq,  (7)  111  pre¬ 
existing  flaws  of  half  length  greater  than  0.04  cm  will  be  activated. 

From  the  flaw  size  distribution  curve  Eq.  (4).  we  calculate  that  there 
are  about  20  flaws  preexisting  in  tho  rock  with  half  lengths  greater  than 
0.04  cm  and  thus  20  cracks  begin  to  propagate.  Since  a  relatively  low 
number  density  of  cracks  is  activated,  crack  growth  will  continue  until 
the  stress  vanishes  (l.e.f  tho  growth  distance  is  not  shortened  by 

coalescence)  and  the  total  fracture  surface  area  produced  is  given  by 

2 

Eq.  (11)  to  be  3.2  cm  .  Based  on  a  qualitative  estimate  of  the  fracture 
damage  in  a  specimen  which  had  experienced  a  similar  stress  history,  this 
estimate  seems  reasonable. 

Finally  it  should  be  pointed  out  that  the  tensile  loading  phase 
experienced  by  a  specimen  which  suffers  fracture  damage  is  not  a  simple 
square  wave.  The  fracture  damage  as  it  develops  modifies  the  tensile 
pulse  continuously,  for  the  individual  cracks  as  they  nucleate  and  grow 
emit  recompression  waves  which  propagate  through  the  specimen  and  unload 
the  material.  Sophisticated  computer  codes  which  calculate  crack  accelera¬ 
tion,  crack  opening,  and  the  stress  histories  in  specimens  undergoing 
fracture  have  been  developed  at  SRI  on  other  projects.  We  expect  to  make 
use  of  these  codes  as  we  further  develop  the  model  for  dynamic  fracture 
in  rock. 

Stage  3:  Crack  coalescence  and  creek  branching — As  radially  expand¬ 
ing  cracks  on  nearby  planes  approach  each  other,  the  magnified  stress 
fields  about  the  crack  tips  overlap  and  the  maximum  stress  direction  is 
altered.  Since  it  is  energetically  most  favorable  for  cracks  to  run 
normal  to  the  direction  of  maximum  tensile  stress,  the  cracks  turn  into 
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each  other,  and  join  up.  Fragment s  are  produced  when  a  sufficient  number 
of  cracks  join  up  to  form  a  continuous  surface  about  an  unfractured 
segment  of  material.  -Figure  23  is  a  photograph  of  a  cross  section  of 
specimen  21  showing  Incipient  In  situ  crack  coalescence  and  fragmentation. 

Another  mechanism  of  fragmentation  is  by  branching  of  Individual 

cracks.  Crack  branching  has  long  been  observed  in  brittle  materials  and 

42 

has  been  extensively  studied  in  glass.  The  phenomenon  is  believed  to 
occur  when  the  crack  speed  reaches  the  maximum  velocity  of  about  c  /3. 
The*-.’  is  evidence  of  branching  in  novaculite  in  Figure  23,  It  is  apparent 
how  this  phenomenon  can  contribute  to  fragmentation  of  the  specimen. 

Stage  4:  Fragmentation — The  duration  of  the  stress  pulse  will 

influence  the  extent  of  fragmentation  produced  by  either  mechanism.  In 

the  present  experiments  stress  durations  of  about  a  microsecond  permitted 

cracks  growing  at  near  theoretical  velocities  in  unfractured  material  to 

reach  a  diameter  of  about  four  millimeters.  This  means  that  at  stress 

2 

levels  sufficient  to  activate  about  five  preexisting  cracks/cm  or  more, 
fragmentation  will  occur. 

The  coalescence  and  fragmentation  phases  of  the  dynamic  fracture 
model  will  receive  considerable  attention  during  the  next  year.  In 
particular  the  fragment  size  distribution  will  be  determined  by  direct 
measurement  of  loost  fragments  and  possibly  in  intact  specimens  by  a 
section  and  transform  method  similar  to  that  used  to  obtain  the  pre- 
existing  defect  size  distribution,  These  results  will  be  correlated  with 
the  applied  stress  and  stress  duration  to  complete  the  predictive  capabil¬ 
ity  for  novaculite. 
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VI  PLANS  FOR  FUTURE  WORK 

The  capability  of  controlling  to  some  extent  the  particle  size 
distribution  produced  by  rock  blasting  would  have  obvious  benefits. 

Two  of  the  most  important  are:  (1)  more  economical  and  efficient  use  of 
charges  to  excavate  a  given  site  and  (2)  production  of  fragment  sizes 
moat  suitable  for  a  given  application  or  for  easy  handling.  Some  ideas 
on  how  to  achieve  such  a  capability  have  been  developed  during  the  course 
of  this  years  work  and  are  discussed  in  this  section. 

Three  factors  are  important  in  determining  the  fragment  size 
distribution  resulting  from  dynamic  loading  of  rock:  the  initial  defect 
structure  in  the  rock,  the  inherent  fracture  resistance  of  crack-free 
material,  and  the  applied  stress  history.  Only  the  last  of  these  can  be 
easily  varied;  quantitative  measures  of  the  first  two  factors  must  be 
obtained.  During  this  first  year  convenient  methods  for  determining 
initial  crack  size  distributions  and  inherent  fracture  resistance  of 
novaculite  were  developed,  and  stress  histories  were  measured. 

The  first  task  in  the  second  year  will  be  to  determine  in  novaculite 
the  fragment  size  distributions  which  result  from  dynamic  loads.  This 
information  will  be  correlated  with  the  initial  defect  distribution,  the 
inherent  fracture  resistance,  and  the  measured  stress  history  to  complete 
the  predictive  capability  for  novaculite.  Following  this,  we  will  begin 
to  test  the  applicability  of  the  model  to  other  rocks. 

The  fragmentation  experiments  on  novaculite  will  be  performed  with 
the  gas  guns  at  various  stress  levels  significantly  above  the  fracture 
strength  to  produce  various  degrees  of  specimen  comminution.  We  will 
develop  techniques  to  recover  100%  of  the  fragmented  specimens,  and 
quantitatively  determine  the  fragment  size  distribution  by  sieve  analyses. 
We  will  also  attempt  to  recover  intact  specimens  in  which  large  numbers 
of  activated  cracks  have  coalesced.  Such  specimens  could  be  analyzed  by 
f ractographic  examination  of  diametrical  sections.  Stress  wave  propagation 
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in  novaculite  will  be  calculated  by  computer  codes,  utilizing  constitutive 
relations  that  include  the  effects  of  stress  relaxation  caused  by  developing 
fracture  damage. 

To  verify  and  generalize  the  dynamic  fracture  model  developed  for 
novaculite,  the  dynamic  fracture  behavior  of  Sioux  quartzite,  Westerly 
granite  and  Tennessee  marble  will  be  studied.  These  rocks  will  be 
characterized  with  respect  to  static  strength,  defect  and  microstructure, 
inherent  fracture  resistance,  and  other  properties  influential  to  the 
fracture  behavior.  Dynamic  tensile  tests  will  be  performed  with  the  gas 
guns  to  obtain  the  dynamic  fracture  strength,  the  inherent  fracture 
resistance  under  high  rate  loading,  and  the  strain  rate  sensitivity. 

Tests  on  both  wet  and  dry  specimens  of  Westerly  granite  will  be  made. 
Ytterbium  stress  gages  will  be  used  to  monitor  stress  histories  of  the 
various  rocks  as  fracture  occurs. 

Fragmentation  experiments  similar  to  those  planned  for  novaculite 
will  be  carried  out  on  Sioux  quartzite  in  the  second  year  and  on  the 
other  rocks  as  time  and  fund  permit.  The  data  will  be  used  to  check 
the  predictions  of  the  model  developed  for  novaculite.  Modifications 
will  be  made  necessary  to  allow  the  model  to  apply  to  rocks  in  general. 
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SUMMARY 

A  simple  homogeneous  rock,  Arkansas  novaculite,  was  chosen  for 

dynamic  tensile  failure  studies.  Experimental  techniques  for  performing 

dynamic  tensile  tests  on  rocks  were  developed  and  used  to  provide  meaningful 

measurements  of  dynamic  tensile  strength,  A  way  to  eliminate  the  occurrence 

of  undesirable  radial  cracking  of  the  specimen  during  testing  was  found; 

specimen  recovery  and  fracture  analysis  methods  were  implemented, 

Novaculite  specimens  were  characterized  with  respect  to  static  and 

dynamic  fracture  strength,  fracture  toughness,  defect  and  microstructure, 

and  other  properties  and  parameters  pertinent  to  the  fracture  behavior. 

A  particularly  important  accomplishment  was  the  quantitative  determination 

of  the  size  distribution  of  inherent  cracklike  defects,  for  this  is 

thought  to  play  a  dominant  role  in  the  comminution  characteristics  of  the 

rock,  and  is  thus  a  necessary  parameter  for  a  dynamic  fracture  model.  The 

quasi  static  tensile  strength  was  determined  by  means  of  expanding  ring 

2 

tests  to  be  440  ±  20  kg/cm  .  The  tensile  fracture  strength  of  novaculite 

2 

was  measured  in  high  loading  rate  gas  gun  experiments  to  be  420  ±  20  kg/cm 

and  found  to  be  independent  of  defect  orientation.  Thus  the  tensile  fracture 

-4  4  -1 

strength  is  strain  rate  insensitive  in  the  range  10  to  10  sec  .  A 
Griff ith-Irwin  type  relation  is  suggested  for  the  failure  criterion  under 
dynamic  loads 


2  -1/2 
a,  =  10  c 
f  max 


where  <r  is  the  fracture  strength  in  kilograms/cm  and  c  is  the  maxi- 
f  max 

mum  crack  radius  in  centimeters . 


The  fracture  mechanism  in  novaculite  under  dynamic  tension  was 
deduced  from  fractographic  observations  to  consist  of  the  following 
sequence  of  events:  (1)  activation  of  preexisting  cracklike  defects, 

(2)  growth  of  individual  cracks  radially  outward,  (3)  coalescence  of 
neighboring  cracks,  and  (4)  fragmentation  of  the  sample.  Based  on  this 
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mechanism  a  dynamic  fracture  model  has  been  partially  developed.  According 
to  this  model,  the  most  influential  parameters  affecting  the  comminution 
characteristics  of  rock  are  the  size  distribution  of  preexisting  cracks, 
the  applied  stress  and  stress  duration,  and  the  Inherent  fracture  resistance 
of  the  rock.  It  is  hoped  that  these  quantities  can  be  related  by  means 
of  the  dynamic  fracture  model  to  predict  resultant  fragment  size  distri¬ 
butions. 

During  this  first  year,  considerable  effort  was  placed  on  determination 
of  the  quantities  associated  with  the  activation  and  growth  stages  of  the 
model.  In  the  second  year,  the  crack  coalescence  and  fragmentation  stages 
will  be  treated,  and  the  dynamic  fracture  model  will  be  completed.  The 
applicability  of  the  model  to  rocks  in  general  will  be  tested  by  perform¬ 
ing  experiments  on  other  rock  types. 
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APPENDIX 


\ 


CALCULATION  OF  THE  PEAK  AXIAL  TENSILE  STRESS  FROM  THE  IMPACT  VELOCITY 
(PLEXIGLASS  IMPACTING  NOVACULITE) 


It  is  assumed  that  both  novaculite  and  plexiglass  are  linearly  elastic 

in  the  stress  region  of  these  experiments  and  that  there  is  no  change  of 

slope  in  the  stress-particle  velocity  curve  going  from  compression  to 

tension.  The  shock  impedances  I  of  novaculite  and  plexiglass  are  then 

s 

calculated  from  the  densities  p  and  sound  speeds  c  acceding  to  the 

o  I 

relation 
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where  the  densities  of  novaculite  and  plexiglass  are  2.63  and  1.18  g/cm  , 

and  the  longitudinal  sound  speeds  at  peak  stress  are  taken  to  be  6.00  and 
3  * 

2.80  x  10  m/sec,  respectively.  The  magnitude  of  k  depends  on  the  units 

-2 

used  for  a  and  u  and  equals  1.02  x  10  when  stress  is  measured  in 
2  X  P 

kg /cm  and  velocity  in  m/sec.  The  shock  impedances  so  calculated  for 

2 

novaculite  and  plexiglass  are  161  and  33.7  (kg/cm  ) (sec/m),  respectively. 
The  peak  axial  tensile  stress  is  given  by 
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since  the  particle  velocity  u  is  related  to  the  impact  velocity  U 

P  proj 

in  the  following  way 
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*  The  density  and  sound  velocity  of  novaculite  were  measured  in  this 
work  as  described  in  Section  III.  The  data  for  plexiglass  are  the 
results  of  Barker  and  Hollenbach,  J.  Appl .  Phys.  4..  4208  (1970). 
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Substituting  the  above  calculated  values  of  shock  impedance  into 
Eq.  (A-3)  we  obtain  the  following  expression, 

'  27'9  Voj  ■  <A-5> 

which  is  used  to  calculate  the  peak  axial  compressive  stress  in  novaculite 
2 

in  kg/cm  from  the  impact  velocity  in  m/sec.  For  the  experiments  with 
a  free  surface,  the  peak  compressivo  stress  should  be  equal  in  magnitude 
to  the  peak  tensile  stress,  assuming  the  tensile  strength  is  not  exceeded. 
If  a  low  impedance  backing  material  is  used,  the  peak  tensile  stress  in 
the  specimen  will  be  equal  in  magnitude  to  the  difference  between  the 
peak  compressive  stress  in  the  specimen  and  that  reached  in  the  backing 
material. 

The  stress  values  calculated,  using  Eq.  (A-5),from  the  measured 
impact  velocities  in  this  work  are  tabulated  in  Tables  II  and  III  and 
are  considered  accurate  to  about  4%,  excluding  effects  of  impact  tilt 
and  material  non-uniformities. 
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